Neither a type-I nor a type-II superconductivity in a two-gap system 
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We show that in two-gap Ginzburg-Landau model there is a wide range of parameters where 
behaviour of a superconductor in external field can not be classified nor as type-I nor as type-II. 
In this regime a superconductor features a nonmonotonic interaction between vortices (which is 
attractive at a certain large length scale and repulsive at a shorter length scale). This results in a 
first order transition into an inhomogeneous state of a mixture of vortex clusters where one of the 
order parameters is suppressed and domains of two-gap Meissner state. 



Interaction of vortices is an interesting question even 
in the simplest one-component Ginzburg-Landau (GL) 
Model . Its character in this model depends on the GL 
parameter k = A/£. In the type-II case (when k > l/\/2) 
the vortices repel each other. Another feature of the 
regime with k > 1/ v2 is that a multiquanta vortex is 
unstable against decay into one-flux quanta vortices. An 
invasion of vortices in a type-II superconductor manifests 
itself as a second order transition at a critical value of ap- 
plied field H c i . In contrast, in the type-I case (k < l/v2) 
vortices attract each other. A system of two of these 
vortices is unstable against a collapse into a single two- 
quanta vortex. In the type-I regime a vortex can not 
form in an applied external field since it would require 
external field which exceeds thermodynamic critical mag- 
netic field. However an isolated type-I vortex is stable. 
The type-I vortices were considered (along with type-II 
vortices) in cosmology Q, where they play important 
role and should be produced during a symmetry break- 
ing phase transition in early universe |3j. The special 
case k = l/y/2 is also very interesting, in particular it 
was shown exactly that the vortices at this value of n do 
not interact p|. 

A general formula for interaction of two parallel vor- 
tices with a separation r is |5( 



U(r) = cl(\,$)K Q 



(i) 



where Kg(r) is the Bessel function, ci,2(A, £) are slow 
functions of the GL parameter, and ci(k = l/y/2) = 
C2(k — 1/V2)- The first term in (JJJ is interaction due to 
particle currents. Indeed this interaction has the range 
of the London penetration A. The second term is the 
attractive interaction due to overlap of cores. 

However experiments |6| show that in certain materi- 
als with k w 1/ vortices attract each other at a large 
scale and repel each other at a shorter distance. The 
nonmonotic interaction of vortices was reproduced in a 
modified one-gap GL functional with additional terms 
in the regime K » 1/ \[2 0, @ • Besides that we should 
mention a nonmonotonic interaction which takes place 
in layered systems: there a van der Waals - type of at- 
traction can be produced by thermal Auctions or disorder 



Below we consider a simplest two-gap model, without 
addition of the extra terms of the type discussed in con- 
text of single-gap model in @ , and without restrict- 
ing discussion to the case k ~ l/v2- We show that 
two-gap GL model allows a novel type of vortices (the 
"double-core" vortices) which have nonmonotonic inter- 
action in a wide range of parameters as their intrinsic 
feature. The effective attraction of these vortices has a 
principally different nature comparing to the effects con- 
sidered in H, @- Moreover, we show that two-gap 
superconductors under certain conditions display even 
more exotic magnetic properties which have no counter- 
parts in single-gap systems. 

The two-gap superconductor 

[HEl (known m par- 
ticle physics as a two-Higgs doublet model 0) in the 
simplest case can be described by the following GL func- 
tional: 
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where * Q = |* Q |e^° and V(|*i, 2 | 2 ) 
Sa=i 2 — M^al 2 + ^I^qI 4 anc > V is a character- 
istic of the interband Josephson coupling strength 



Recently two-gap superconductivity was observed in in 
MgB2 13] and in NbSe2 [13 ■ Two-gap models appear 
also in the theoretical studies of liquid metallic hydrogen, 
which should allow superconductivity of both electronic 
and protonic Cooper pairs |l5j|. while in liquid metal- 
lic deuterium a deuteron superfluidity may coexist with 
electronic superconductivity [l^. Besides that in super- 
conductors there could coexist two types of pairing (e.g. 
s- and p-wave condensates) . 

One can rearrange variables in (J2J) as follows [lfiL ItR 



2 



0: 



2 |^i| 2 |1>2| 2 , h 2 

mi r?i2 mi rn-2 

Jl*i| 2 | $2 |2 2 Y 

< Vcpi + V02 A } 

Ami 4m2 c 

H 



+ 2»j|*i*2| cos(^ - <h) + V(|*i, 2 | 2 ) (3) 

Let us consider first the case rj = 0. The first term fea- 
tures a gauge invariant phase difference 7„ = <f>\ — <p2- 
Since no coupling to gauge field enters this term, then 
any vortex with nontrivial winding number of 7„ has di- 
vergent energy per unit length |l7|. In this paper we 
consider only vortices with finite energy per unit length. 
These are the vortices with phase windings satisfying the 
condition A^ n = A(<pi — fc) = where Aj n denotes the 
gain in the phase "f n when we go around a vortex core. In 
the case of nonzero 77, one can observe that, when rj > 
the free energy J3J is minimized by <f>\ — 4>2 = 7r, while 
when 77 < we have <j)\ — 4>2 = 0. Thus when one con- 
siders topological defects with A^f n = the fourth term 
in is phase-independent and can be absorbed in the 
potential term. In this paper we want to discuss mag- 
netic properties of two gap system in the simplest and 
most transparent case and for this reason we put 77 = 0. 
As discussed above the low-temperature magnetic prop- 
erties of the cases of zero and nonzero 77 are essentially 
the same since in both cases only the vortices character- 
ized by A-jn = have finite energy. An extension to a 
nonzero rj as well as to an addition of terms proportional 
to I 'I'll 2 1 ^2 1 2 is trivial. So for the case when A-y n = 0, 
the remaining nontrivial terms in GL functional are: 




where 7 C = (</>i + < t> 2 ) The equation for supercurrent for 
is (see e.g. El [Hill) 
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As it was shown in |17| . in this model, the vortices char- 
acterized by the following phase gains: A"f n = 0; A'jc = 
Airk carry /c-quanta of magnetic flux. 

The model (@J possesses three characteristic length 
scales: two coherence lengths £ Q = hf \J<±vn a b a deter- 
mined by the effective potential and the magnetic field 
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Let also introduce a notation X a = c^/m a j (v87re|\& Q |) 
Below we also use notations for thermonynamic criti- 
cal magnetic fields for individual condensates £f c t( Q ) = 
$>o / (2y/2ir{; a \ a ) . In discussion when the individual con- 
densate *i is of type-II we denote the first and the sec- 
ond critical magnetic of isolated \&i as follows: H c im — 
$ /(47rA 2 )[log(Ai/a) + 0.08] and H c2(1) = $ /(2< 2 ). 

Apparently in the cases when £1 « £2 >> A and when 
£1 ~ £2 << A, magnetic properties of the two-gap super- 
conductor parallel those of single-gap type-I and type-II 
superconductors correspondingly. However in a general 
case, when one of the condensates is of type I and another 
is of type II, a two-gap system in exernal field has much 
richer phase diagram and allows many types of different 
transitions in external field. 

In the case when £1 << Ai and £2 >> A2, and when 
thermodynamic critical magnetic field H c t{2) for the type- 
I condensate 'J 2 is much higher than i? C 2(i) for the type-II 
condensate *i, the system undergoes a first order tran- 
sition from a normal state immediately into a two-gap 
superconducting state. This is because, in this case, 
at fields higher than H ct ( 2 ) the system does not allow 
nontrivial solution of the linearized GL equation for *i, 
while, when the field is lowered below ff c t( 2 )j there ap- 
pears a transition immediately into two-gap supercon- 
ducting state. In this state the magnetic field is screened 
in bulk of a sample largely due to surface current \&2 and 
nothing can preclude the appearance of the second con- 
densate *i in the bulk of the sample even if the applied 
field is much larger than ff C 2(i)- 

Consider now the case when Ai/£i > l/V^',M/^2 > 
1/V2; a « 6; A - (I/A 2 + l/A 2 )- 1 ^ < 6/V2. In this 
case we have a two-gap system made up of two type-II 
condensates, however, somewhat counterintutively, the 
magnetic properties of this system are principally differ- 
ent from a type-II superconductor. One can observe that 
the resulting penetration length in this system can be 
smaller than one of the coherence lengths A < £ 2 /\/2, 
since both condensates contribute to A, so a vortex in 
this system will consist of an inner core of characteristic 
size £1, a flux tube of size A, and an extended "outer" 
core of a characteristic size of £2 larger than A ^3- ^ 
means that, albeit at length scales smaller than A the 
interaction of two vortices is repulsive, at the same time 
the extended core gives to such a vortex an effective at- 
traction U oc Kq[t /^/2^2) for r > A in a range of \/2£. 
Thus in contrast to an ordinary type-II superconductor, 
there exists a range of parameters where the phase tran- 
sition at H c \ in this system is of the first order. This 
is because due to attractive interaction the energy of a 
lattice made up of n vortices (with the lattice spacing de- 
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termined by the minimum of the nonmonotonic interac- 
tion potential) has lower energy than n isolated vortices. 
Correspondingly the magnetization of such a supercon- 
ductor should change abruptly by a formation of a vortex 
lattice, instead of a gradual invasion of isolated vortices. 
Also for the same reason in a decreasing applied field 
the vortex lattice should survive even in fields lower than 
the critical field where a formation of an isolated vor- 
tex becomes energetically favored. Another consequence 
of the effective attraction between vortices is that at a 
given low vortex density, when average distance between 
the vortices is larger than the distance determined by 
the minimum of the effective potential, the vortices will 
spontaneously form lattice clusters in between areas of 
vortex-free Meissner state. 

Let us now turn to the case when Ai/£i > 
< l/\/2;C2 > A. Let us consider a vor- 
tex in this case. For a two-gap superconductor, vor- 
tex energy consists of energies of cores, kinetic energy 
of Meissner current and magnetic field energy. The 
magnetic field energy and the kinetic energy of the 
screening current is given by the standard expression: 
F m = (1/8tt) J d 3 xH 2 + (l/8vr) J d^A^curlH) 2 . Here 
we note that in the present case | \&2 (^) | 2 varies slowly 
over the London penetration length A. Magnetic field is 
screened at a distance from the core which is smaller than 
£2- That means that only depleted density of Cooper 
pairs of condensate $2 participates in the screening of 
magnetic field. Thus one can not use the London pene- 
tration length A but should introduce an effective pene- 
tration length A e// = [l/Ai + i/A!(x)]~ 1/2 where A 2 (x) = 
c v /m2"/(v / 87re|* 2 (x)|) > A 2 = c v ^ 2 "/(v / 87re|* 2 |), where 
1^/2 1 is the average value of the order parameter mod- 
ulus. In the case when Ai << £2, we have A e // « Ai 
which corresponds to the situation when magnetic field 
is screened mostly by condensate ^1. On the other hand 
in the case when £ 2 >> Ai >> A 2 >> £1 the magnetic 
field is screened at a chracteristic length scale of Pippard 
penetration length A = Af « (A^) 1 ^ 3 of the conden- 
sate \E , 2-This expression is valid when A^ < Ai. In the 
above expression for F m we cut off integrals at the dis- 
tance £1 from the center of the core in order to obtain 
an estimate of the vortex energy which has a logarith- 
mic accuracy. With it the energy per unit length of a 
onc-flux-quanta vortex is 

^(j^r-) 2 ^^ + v cl + v c2 (7) 

\4mX eff J £1 

Where V ca arc the energies of the cores per unit 
length which are of order of magnitude of [core size] x 
[condensation energy] . The estimate of the core energy 
can also be expressed as: 

v jf g^ _ 1 /gp\ 2 1 

Vca ~ 8tt ~ Sir c 2 m a ~ 4 {An J A 2 { ' 



Apparently V c i + V c2 « ($ /8ir\) 2 . 

A straightforward calculation of the field charac- 
terizing when it is becoming energetically preferred to let 
into a superconductor a single vortex gives: 
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The transition into mixed state should be of the first or- 
der since the extended core at length scale £2 > A e // 
gives rise to attractive interaction of vortices so, a lattice 
of vortices with a spacing determined by the minimum 
of the interaction potential is preferred over a system of 
widely separated vortices. That is, the energy of a system 
of n vortices is minimized when vortices spontaneously 
form a lattice cluster with overlapping outer cores. We 
should observe that in the present situation is larger 
than thermodynamic critical magnetic field of the type-I 
condensate which is H ct ( 2 ) = $o/(2v / 2ttC2A 2 ). We stress 
that this however does not mean that at the field the 
condensate ^ 2 is completely depleted due to overlapping 
of outer cores. It is because when the applied field is 
close to H^, then the field is mostly screened by a su- 
percurrent of the condensate which circulates along 
sample's edge. Thus the vortex system is dilute, or more 
precisely the intervortex distance is only determined by 
the effective attraction. So in contrast to original type-I 
and type-II behavior, the two-gap superconductor in this 
regime displays a first order transition into an inhomo- 
geneous state consisting of clusters of vortices, where the 
order parameter ^ 2 is depleted due to overlaps of outer 
cores. So the superconductivity in these vortex lattice 
"droplets" is dominated by the order parameter ^1 and 
is essentially a one-gap superconductivity. Since the vor- 
tex density depends on applied field while the intervor- 
tex distance is determined by the nonmonotonic interac- 
tion potential, so, besides these clusters of vortices, there 
should be present domains of two-gap superconductivity 
in the vortexless Meissner state. 

One of the possible direction where this work can be 
extended is the interesting question for a numerical study 
if there is a range of parameters when a triangular lat- 
tice is not the most energetically preferred. Indeed, for a 
given density of vortices, at a certain applied field, spac- 
ing of a square lattice may be closer to length scale given 
by minimum of vortex interaction potential, than spac- 
ing of a triangular lattice. However, in a triangular lat- 
tice, albeit the intervortex distance is larger, there are 
six instead of four nearest neighbors. Apparently in suf- 
ficiently strong magnetic fields the effective attraction at 
large length scale is irrelevant and vortices form a stan- 
dard triangular lattice. So this system may allow tran- 
sitions between lattices with different symmetries, when 
applyed magnetic field is varied. 

We also would like to note that in the recent experi- 
ment on MgB 2 |20j (which is the strictly type-II two-gap 
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superconductor) indicates a large disparity in coherence 
lengths. The coherence length in 7r-band was found being 
^50nm which is much larger than an estimate which one 
would obtain from a standard GL formula for a one-gap 
superconductor using the H C 2 value of MgBi. Such a 
behavior is completely in agreement with the model for 
the composite two-gap vortex Pjfl ■ 

In conclusion we discussed magnetic properties of a 
simplest form of a two-gap GL model. We have shown 
that the model exhibits behavior which has no counter- 
part in a standard one-gap superconductors. The dis- 
cussed features do not allow to classify the two-band su- 
perconductor in the considered regimes as type-II and 
indeed nor as type-I but it should legitimately be distin- 
guished in a separate class. 



Note added A more detailed and quantitative study 
of the question which we discuss in this article can be 
found in E. Babaev, J.M. Speight cond-mat/0411681 



[1] A. A. Abrikosov Sov. Phys. JETP 5 1174 (1957). This 
concept was introduced in high energy physics in 
H. B. Nielsen and P. Olesen, Nucl. Phys. B 61, 45 (1973). 

[2] M. B. Hindmarsh and T. W. Kibble, Rept. Prog. Phys. 
58, 477 (1995) 

[3] T. W. Kibble, J. Phys. A 9, 1387 (1976). 

[4] E. B. Bogomolny, Sov. J. Nucl. Phys. 24, 449 (1976) L. 
Jacobs, C. Rebbi Phys. Rev. B 19 4486 (1979). 

[5] L. Kramer Phys. Rev.B 3, 3821 (1971). 

[6] R. P. Hiibener Magnetic Flux Structures in Superconduc- 



tors Springer, Berlin (1979); U. Krageloh Phys. Status. 

Solidi 42, 559 (1970). 
[7] A. E. Jacobs, Phys. Rev. B 4, 3022 (1971), A. E. Jacobs, 

Phys. Rev. B 4, 3029 (1971); A. Hubert Phys. Stat. Sol. 

(b) 53 147 (1972); E.H. Brandt ibid 77105 (1976); I. 

Luk'yanchuk Phys. Rev. B63, 174504 (2001). 
[8] F. Mohamed, I. Luk'yanchuk, M. Troyer, G. Blatter 

Phys. Rev. B 65 224504 (2002). 
[9] G. Blatter and V.B. Geshkenbein, Phys. Rev. Lett. 77, 

4958 (1996). S. Mukherji and T. Nattermann, ibid 79, 

139 (1997). 

[10] H. Suhl, B.T. Matthias, L.R. Walker, Phys.Rev.Lett.3, 
552 (1959) 

[11] A.J. Leggett, Prog, of Theor.Phys. 36, 901 (1966) 

[12] T.D. Lee, Phys.Rev. D 8 1226, (1973) 

[13] see e.g. F. Bouquet et. al. Phys. Rev. Lett. 87 047001 
(2001) Amy Y. Liu et. al, Phys. Rev. Lett. 87 087005 
(2001) P. Szabo et. al. Phys. Rev. Lett. 87, 137005 (2001) 

[14] T. Yokoya et. al. Science 294 2518 (2001) 

[15] K. Moulopoulos and N. W. Ashcroft Phys. Rev. B 59 
12309 (1999); J. Oliva and N. W. Ashcroft Phys. Rev. B 
30 1326 (1984); N. W. Ashcroft J. Phys. A129 12 (2000). 

[16] E. Babaev, L. D. Faddeev and A. J. Niemi, Phys. Rev. B 
65, 100512 (2002) Similar consideration were found be- 
ing relevant for spin-triplet superconductors E. Babaev, 
Phys. Rev. Lett. 88, 177002 (2002) 

[17] E. Babae v, Phys. Rev. Lett. 8 9, 067001 (2002) 

[18] E.Babaev cond-mat/0201547 

[19] In fact the actual penetration length will be larger than 
the London penetration length A since it will receive a 
correction due to depletion of |SJ/2 1 at a distance smaller 
than £2 from the core. The effective penetration length 
A e // in this vortex, depending on the GL parameters, 
will have the range A < A e // < Ai < £2 

[20] M. R. Eskildsen et. al. Phys. Rev. Lett. 89 187003 (2002) 



